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ABSTRACT 



Context. Hydrides of the most abundant heavier elements are fundamental molecules in cosmic chemistry. Some of them trace gas 
irradiated by UV or X-rays. 

Aims. We explore the abundances of major hydrides in W3 IRS5, a prototypical region of high-mass star formation. 

C/2 Methods. W3 IRS5 was observed by HIFI on the Herschel Space Observatory with deep integration (~ 2500 s) in 8 spectral regions. 

^ C^^ Results. The target lines including CH, NH, H3O+, and the new molecules SH+, H2O+, and 0H+ are detected. The H2O+ and 

" 0H+ J = 1 — lines are found mostly in absorption, but also appear to exhibit weak emission (P-Cyg-like). Emission requires 

, I high density, thus originates most likely near the protostar. This is corroborated by the absence of line shifts relative to the young 

►^ , stellar object (YSO). In addition, H2O+ and 0H+ also contain strong absorption components at a velocity shifted relative to W3 

JI^ ■ IRS5, which are attributed to foreground clouds. 

1^^ ■ Conclusions. The molecular column densities derived from observations correlate well with the predictions of a model that assumes 
1— ' the main emission region is in outflow walls, heated and irradiated by protostellar UV radiation. 

CO ' Key words. Stars: formation - stars: high mass - ISM: molecules - ISM: individual objects: W3 IRS5 - Line: identification - 
►^ ■ Ultraviolet: ISM 

O. 

o. 

^— H ■ 1. Introduction (FUV) or X-rays - interact with the molecular gas and heat 

IL* ^ . „ 1 T 1 . 1 . • 1 1 T it, hydrides and particularly their ions are greatly enhanced 

^> In interstellar clouds, chemical reactions with hydrogen .^ abundance (e.g., Hollenbach & Tielens 1999). Ionized 

K^ molecules lead to an elementary class of molecules that ^^^.^^^ ^^^ chemically even more active and can substan- 

^ represent key species m the chemical evolution to larger ^.^^ ^^.^^ chemical evolution. Many of the above hydrides 

_ C^ _ molecules. These fundamental mokcules, known as hy- j^^^^ ^^^^ observed from the ground in absorption in diffuse 

iVr+''Z^i c,?-H 't?^^ ' h £' ""^u ^^^^"^ 3°^^' ^^ i interstellar clouds (e.g. Swings & Rosenfeld 1937; Menten 

CH+, NH+, SH+, H2O+, and H3O+. The combmation of ^^ ^j 2OIO; Wyrowski et al. 2010a). 

hydrogen atoms with a heavier atom causes large dipole mo- ^^ ^ i i • t • i . r • 

, -11 ,,. ,, X-11--1-- • Here we locus on hydrides m dense star-iormmg re- 

ments and large rotation constants, particularly m diatomic . ,. ,^, ^^^ . n . if r 

u J -J rj^u- -J 1 4- ^u --L ^- 1 1 r-k 1 gions surrounding protostars. We report nrst results oi 

hydrides. 1 his widely separates the excitation levels. Only ^, .t^ t . • i-f- .• , .,.,.. ,, tt ^ ^ 

low- J lines are excited at temperatures relevant to star and f^« Radiation Diagnostics project withm the Hersche 

1 , r -L- rr^u 1- 1, u u ui Kcy program Water iu Star- tormiug rcgious With iier5c/ie/ 

planet formation. These lines have now become observable .^^^^^^^^^ ^. , , ^ , o^^r^x rr^i .t-* t .• 

•-H. ^u u h 1 Q niu ^ ^o-iu 4-^ ^ 1 onin\ (WISH, van DishoecK et al. 2010). Ihe Radiation 
with the iier5c/ie/ Space Observatory (Filbratt et al. 2010). \. .. , . . . . i- i t 
7v/r uj-j u u-u 4-- \- ■ 4-u ■ Diagnostics subprogram aims to trace high-energy radia- 
Many hydrides have a high activation energy m their ,. ,. i . i . i ^ , ^ i • i i . i • • i 
r -L- -lu Tj -fu- u u -^ f jj^T tion through the most abundant hydrides and their ions and 
formation paths. However, it high-energy photons - tar UV 1..1..1 1 -i . ir . t^ • r 
relate them to the chemical network of water. Preparing for 

* Herschel is an ESA space observatory with science instru- this project, Stauber et al. (2005, 2006) showed that hy- 
ments provided by a European-led Principal Investigator con- drides may be enhanced in star-forming regions affected by 
sortia and with important participation from NASA far UV and/or X-ray irradiation. In a second step, optimal 
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lines were selected using catalogued data in the available 
databases or computed from measured molecular constants 
(Bruderer 2006). 

Spherically symmetric models including FUV and X- 
rays have found that most hydrides would be too weak for 
Herschel detection. However, evidence of extended high- 
energy irradiation (Doty et al. 2004) and in particular the 
observation of CO^, with an abundance four orders of mag- 
nitude larger than predicted, in W3 IRS5 and other YSOs 
with the JCMT (Stauber et al. 2007) raised the expecta- 
tion that the effects of irradiation may be more dramatic 
in asymmetric reality. Bruderer et al. (2009) modeled and 
interpreted the C0+ anomaly as FUV radiation originating 
in the YSOs and irradiating the walls of outflow cavities. 

Based on these predictions, the 'Radiation Diagnostics' 
observations started with deep integrations of a large num- 
ber of hydrides that was to be followed by a survey of a few 
species in many sources of various ages and masses. Here 
we present exploratory observations towards W3 IRS5, a 
nearby region (1.83 kpc, Imai et al. 2000) of high-mass star 
formation, moving at - 38.4 km s~^ in the local standard 
of rest (LSR). W3 IRS5 resembles the Trapezium cluster in 
Orion but is considerably younger. At least six radio sources 
are in the Herschel beam. They represent 'hyper compact' 
H H regions produced by high-mass YSOs, at least two of 
which are O stars (van der Tak, Tuthill, & Danchi 2005; 
Rodon et al. 2008). In this paper, we report the observa- 
tions of major hydrides towards W3 IRS5 apart from H2O, 
CH+, and OH, which are studied in other projects of WISH 
(Chavarria et al. 2010; Wampfler et al. in prep.). 



2. Observations 

The Heterodyne Instrument for the Far Infrared (HIFI, de 
Graauw et al. 2010) on iJer^c/ie/ observed W3 IRS5 between 
1 and 8 March 2010 in the Science Demonstration Phase 
in eight 4 GHz frequency bands for about 2500 s each. One 
of them includes [C II] at 1900.5369 GHz. We used the 
Wide Band Spectrometer which has a spectral resolution 
of 1.1 MHz, and HIPE 3.0 for pipeline and data analysis. 
The data were taken by double beam switching (DBS), the 
high-frequency (HEB) bands in fast DBS mode. The off- 
source position was at a distance of 3 arcmin in the NE 
and SW without remarkable IR sources. 

The current accuracy of the velocity calibration is es- 
timated to be better than 2 km s~^. The antenna tem- 
perature was converted to main beam temperature, using 
pre- flight antenna efficiencies. After visual inspection and 
defringing, the V polarization was shifted linearly in flux to 
match the H polarization. The two polarizations were then 
added. A second observation of equal length was made using 
a local oscillator frequency shifted by 10 km s~^. The two 
data sets were plotted in both upper and lower sideband 
presentation. If a line matched velocities in one sideband 
and was double with 20 km s~^ separation in the other, 
the frequency of the former was assumed. All lines of inter- 
est could be attributed to a sideband without ambiguity. 
The continuum was divided by 2 for double sideband ob- 
servations, assuming that it is the same in both sidebands, 
of equal sensitivity. 



3. Results 

Most lines (except H3O+) are split by fine or hyperfine 
interaction as indicated in Fig. 1. Table 1 lists the ob- 
served lines (strongest only for multiples) and summarizes 
the quantitative observational results. Molecules here de- 
tected for the first time in star-forming regions include 
H2O+, 0H+, and SH+. Having the most prominent line 
near the H2O para ground-state line at 1113.3 GHz, H2O+ 
is serendipitously detected in many Herschel observations. 
H2O+ and 0H+ are detected in absorption by the inter- 
stellar medium (Bruderer et al. 2010b; Ceccarelh et al. 
2010; Gerin et al. 2010; Falgarone et al. 2010; Neufeld et 
al. (2010); Ossenkopf et al. 2010; Schilke et al. 2010), but 
also near the systemic velocity of other high-mass YSOs 
(Wyrowskiet al. 2010b). 

Several lines are found in absorption, indicated in Table 
1 by negative peak values relative to the continuum. All 
lines predominantly in absorption (NH, 0H+, H2O+) orig- 
inate in molecules in the ground state with a J = 1 level 
energy exceeding 47 K. Lines of CH and SH+, observed in 
emission, are transitions from a J = 1 level of energy less 
than 26 K. We do not detect SH and NH+ above the 5<j 
limit. 

The observed line shifts are within ±2 km s~^ of the 
YSO, thus within the accuracy of the frequency calibra- 
tion, spectral resolution, and molecular data. The only de- 
viations in Table 1 are 0H+ and H2O+. Their strongest 
peaks, both in absorption, are shifted by 34 - 39 km s~^ 
relative to the YSO, similar to an absorption feature in [C 

II]. 

Figure 2 illustrates a possible explanation of the ex- 
ceptional line shifts of 0H+ and H2O+. In both cases, 
a theoretically strong hyperfine structure line corresponds 
to a smaller absorption dip within less than 5 km s~^ of 
the YSO. For 0H+, the low- frequency line (a doublet) is 
blended with the high-frequency line of another, red-shifted 
0H+ component. We thus interpret the spectra of 0H+ and 
H2O+ as the superimposition of a component originating 
in the star-forming region and another component at lower 
frequency related to an outflow or the foreground interstel- 
lar medium near zero velocity in the LSR. 

The line widths in Table 1 are between 4 and 7 km 
s~^, measured at the component with the systemic velocity 
of the YSO. Two lines attract attention: {i) CH has an 
absorption feature near the YSO velocity. This absorption 
reduces the peak flux and widens the line profile at half 
power, (n) The line width of NH is extremely narrow in all 
hyperfine structure components. 

The lines are generally symmetric. A remarkable excep- 
tion is 0H+, showing blue tails. We note that the shifted 
lines of 0H+ (unblended line) and H2O+ are 17.7 and 11.1 
km s~^ wide, respectively, indicating that they have a dif- 
ferent origin from the component at zero systemic velocity. 

Several lines in Fig. 1 show slightly blue-shifted absorp- 
tion and red-shifted emission, thus a P-Cyg-like behavior. 
Figure 2 compares them with the [C II] line. CH (strongest 
component) exhibits similar absorption to [C II], their peak 
absorptions being at -1.6 km s~^ and -1.7 km s~^, respec- 
tively, relative to the YSO motion. 0H+ and H2O+ differ, 
having absorption peaks at -2.5 and -4.8 km s~^, respec- 
tively, and emission peaks at 4.2 and 3.5 km s~^. 

Column densities are given in Table 1 for the upper 
energy level of transitions in emission, and for the lower 
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Table 1. Frequency, upper level energy, and Einstein coefficient of molecules and lines observed by Herschel/HIFI 
towards W3 IRS 5. The numbers in parentheses give the decimal power. Negative peak fluxes signify background minus 
line temperature of lines in absorption. Line widths refer to the FWHP value of the most intense line peak or absorption 
and its line shift to the systemic velocity of -38.4 km s~^. Non-detected peak fluxes are 5a upper limits at 1 km s~^ 
resolution, and non-detected line fluxes assume a 5 km s~^ line width. The column densities refer to the upper (emission) 
or lower (absorption) state. 
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Fig. 1. Spectral lines observed towards W3 IRS5 in main beam temperature. The positions of the theoretical 
fine/ hyp erfine structure lines are shown shifted by -38.4 km s"-*^, the systemic motion of the YSO. The length of 
the vertical bars indicates the theoretical intensities in arbitrary scale. The frequency binning is 1 km s~^. 



energy level of lines in absorption. The fine /hyp erfine struc- 
ture lines were summed over all components using their sta- 
tistical weight (see Bruderer et al. 2010b for details). The 
observed parameters of H3O+ at 1031 GHz and 1069 GHz 
are extracted from blends (see Fig. 1) by fitting Gaussians. 

A rotational diagram can be constructed for H3O+ (see 
Appendix). The data can be fitted with a straight line, sug- 
gesting that (i) the ortho-to-para ratio is compatible with 



unity and (ii) a unique rotational temperature of 239 K (er- 
ror range 205 - 285 K) exists. It indicates that the H3O+ 
emission originates, on weighted average, from regions with 
equal or even hotter kinetic temperature, or, if pumped, ra- 
diation temperature. 
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Fig. 2. Light hydrides with P Cygni type absorption fea- 
tures compared with [C II] . The flux scales are reduced and 
shifted for better visibility. The velocity is relative to the 
systemic velocity of W3 IRS5, shown by a vertical dashed 
line. The horizontal dashed lines indicate the continuum 
level. 



4. Discussion and conclusions 

The observed shifts of the line peaks relative to the systemic 
velocity of the YSO are small and suggest that the lines 
originate in the star-forming region, not in the foreground 
interstellar medium. This may not be the case, however, for 
the components 0H+ and 112© + shifted by a larger amount 
to roughly to the velocity Vlsr = 0- 

Additional evidence of the origin comes from the lines in 
emission, which are produced at critical densities (if known) 
of order > lO^cm"^. Even lines predominantly in absorp- 
tion, such as 011+ and possibly II2O+, have a small emis- 
sion feature in the red wing, thus a P-Cyg-like profile, near 
the systemic velocity. This is the first report of Galactic 
011+ and H2O+ in emission. Figure 2 shows the similar- 
ity of the profiles of 0H+ and II2O+. The absorption in 
CH may be caused by self-absorption, but is similar to that 
in [C II], which is found in emission at the off position. 
Absorption lines in star-forming regions need detailed mod- 
eling and radiation transfer calculations beyond the scope 
of this letter. 



The ground-state column densities of the component 
moving with the YSO are 8.1 x 10^^ cm~^ for H2O+, and 
9.7 X 10^2 cm-2 for 0H+. The values for the red-shifted 
components are 3.8 x 10^^ cm~^ and 6.1 x 10^^ cm~^, re- 
spectively. Both components yield larger OH+/II2O+ ratios 
than the other observations reported in this volume. 

The measured line widths are generally small (< 7 km 
s~^) and show no anisotropics. We thus find no evidence of 
shocks except possibly in 011+ , which needs to be studied 
in combination with shock tracers. 

In Fig. 3, column densities are displayed, derived from 
integrated line fluxes neglecting re-emission or reabsorption 
of the final state. The derivation is based on Table 1, except 
for 011+ and H2O+ where only the unshifted component 
attributed to the YSO is used. For II30+, the value of the 
J = 4 level as well as the one from the rotational diagram, 
integrated over all levels, are shown. The column densities 
are compared with abundances predicted by Bruderer et al. 
(2010a) in a two-dimensional 'standard' YSO model used 
here as a template, assuming UV and X-ray irradiation by 
a central high-mass YSO. It enhances the abundance of 
diatomic hydrides in the outflow walls by many orders of 
magnitude, such that the beam- averaged abundance is sig- 
nificantly changed. It is averaged over a radius of 20 000 
AU, or 10.9'' at the distance of W3 IRS5. 

The similarity of observations and model abundances in 
Fig. 3 support the scenario of hydride enrichment in outflow 
walls heated and irradiated by protostellar far UV. 
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Fig. 3. Comparison of low- J level column density with pre- 
dicted beam-averaged abundance relative to H2 from a 2D 
model by Bruderer et al. (2010a). 
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Fig. A.l. Rotational diagram of H3O+. Numbers indicate 
the frequencies in GHz of the observed hues. 



Appendix A: 

The rotational diagram in Fig. Al is complemented with 
ground-based data for the frequency range 300 - 400 
GHz observed at the CSO with comparable beam size by 
Phillips, van Dishoeck & Keene (1992). The non-detection 
at 307 MHz is surprising, but possibly an effect of optical 
depth. The data in Fig. A.l, except 307 GHz, are well fitted 
by a single rotational temperature of 239 K, suggesting that 
the observed levels are populated according to an exponen- 
tial distribution. The derived temperature and column den- 
sity infer an optical depth of r < 0.1 for all lines except at 
307 GHz. The fitted line (dashed) corresponds to a column 
density of 8.5 (±2) x 10^^ cm~^, consistent with the value de- 
rived by Phillips, van Dishoeck & Keene (1992). This leads 
to a beam averaged H3O+ abundance of 4.2(±1) x 10~^^ 
relative to H, to be compared with the theoretical value of 
4 xlO~^^ reported by Bruderer et al. (2010a). 



